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ABSTRACT 

This report presents the results of the study of surface elevation 
deviations in the Gulf Stream region off the eastern coast of the United 
States between Wallops Island, Virginia and Miami, Florida. The main causes 
of surface elevation deviations are geoid perturbations due to the continental 
shelf and the gcostrophic adjustment of the density field due to the Gulf 
Stream. Quantitative surface elevation profiles were calculated based on 
geophysical measurements of gravity anomalies and hydrographic data. The 
results are presented graphically along with contemporaneous weather data. 
Comparisons are made between the profiles based on hydrographic data and a 
mean theoretical model. Tlie agreement is generally good. The theory of 
geostrophic flows Includirg some classical Gulf Stream models is also presented 
briefly. 
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1, INTRODUCTION AND SUMMARY OF RESULTS 

1.1 INTRODUCTION 

The objectives of the Ct.lf Stream model study are (a) to find ouc the 
locations and magnitudes of surface elevation deviations due to water movement 
and other causes, (b) to use such results to Interpret the satellite altimeter 
data, and finally, (c) to calculate water movement based solely on satellite 
altimeter measurements. Due to the limiting scope of the present phase of Che 
study, this report presents the results of a study of surface elevation 
deviation due to water movement and gravity anomalies along the Gulf Stream 
between Miami, Florida, and Wallops Island, Virginia. 

In Chapter 2, a brief discussion of the relation between surface elevation 
deviation and the water movement Is presented. This Includes theory of 
geostrophic flow, numerical procedure used to calculate Che surface elevation 
based on hydros’’^ Hara. classical Gulf Stream models, and a mean numerical 
model of the Gulf Stream using all Che existing hydrographic data In the North 
Atlantic Ocean. 

Chapter 3 presents all of the surface elevation profiles calculated by 
the method discussed In Chapter 2 and based on data published by Che National 
Oceanographic Data Center, NOAA, covering all of the ship bound surveys from 
19l2 to 1972. Together with the profiles, wind data are presented whenever 
available. The profiles are also compared with the results of the mean 
numerical model of Che Gulf Scream. The agreement Is generally, and In a mean 
sense, good. 

Chapter 4 deals with the calculation of gcold perturbation at Che 
continental shelf. Due to the lack of data, an approximation method has been 
used. The procedures of approximation, calculation and the results are 
presented. 
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1.2 CONCLUSIONS 

The results of this study Indicate that: 

(a) A mean Gulf Stream model can be constructed but such a model can only 
be used as an estimate of the location and the strength of the Gulf 
Stream. Individual cases may deviate substantially from such a model 
by perhaps lOOZ In strength as a result of seasonal fluctuations and 
also due to local meteorological conditions. 

(b) The mean location of the Gulf Stream south of Cape Hatteras generally follows 
the continental shelf break where large geold perturbations exist due to 

the continental rise. The change of the geold elevation varies from one 
location to another, but. In general, the magnitude Is larger than that 
produced by the Gulf Stream. However, the location and the magnitude of 
the geold perturbation are fixed In space and time. Such features can be 
easily filtered out o£ ubsei. vatlwuul data. 

(c) Due to the high mobility of the Gulf Stream, shlpbound surveys alone are 
definitely Insufficient to determine the nature of the motion. Remote 
sensing methods have to be used, not in lieu of, but rather in complement 
to the standard hydrographic method. 


. r 



1 


1 


} 


3 


2. THEORY AND CALOJUTION OF DYN/\MIC HEIGHTS AND GULF STREAM MODELS 
2,1 INTRODUCTION 

The distribution of temperature, salt, density and pressure within and 
aotion of any body or part therein of water is approximately governed by a 
system of seven non-linear partial differential equations, l.e., 


^ f - -II- - 3f 
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(2.4) 


(2.5) 

( 2 . 6 ) 

(2.7) 


where u is the velocity, A the earth's rotational vector, P the pressure, 
g the gravitation acceleration, f the density, the Kronecker delta, S 

the salinity, T the temperature, >*■ the molecular viscosity, fc the cyclic 
tensor, D the eddy diffusivlty, K the eddy conductivity, ol the coefficient 
of thermal expansion, ^„the coefficient of volumetric expansion and the 
subscripts i, J, k arc each equal to 1, 2 and 3. These equations describe the 
behavior of a taroclinlc fluid in an arbitrarily shaped basin on a spherical 
earth with boundary, bottom, side and surface, conditions which are consistent 
with the physical effects and subsequent constraints of the surrounding environ- 
ment. 
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The complete system of equations (2.1 through 2.7) Is impossible to solve 
given the present state of mathematical art. Under certain conditions however, 
simplifications based on appropriate physical considerations can be made to 
reduce the equations to a more manageable, tractible form. In the case of 
major ocean current systems, the motions arc rather steady and to a lower order 
approximation, it can be assumed that the ratio of the non-linear acceleration 
to Che rotational force, l.e., the vertical component of relative vorticity, is 
rather small. Since molecular friction and vertical accelerations can always 
be neglected in the steady state and furthermore because we may assume that 
turbulent stresses change slowly over distances which are short compared to the 
characteristic length of the flow, then \e are permitted to simulate a particular 
class of motions called "geostrophlc currents." As long as the current system 
under investigation is not in a lateral frictional boundary layer, where sharp 
horizontal gradients of velocity can be generated or in either one or two 
vertical frictional boundary layers which Include strong vertical velocity 
gradients, the geostrophlc assumption can yield useful approximations to the 
actual current system. In such cases, the governing momentum equations (the 
reduced Navler-Stokes' system) become: 
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where the subscript H denotes horizontal components. 

We can integrate (2.9) in z and combine with (2.8) to get 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 
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where Is the deviation of the free surface from mean sea level. Obvious 
factors Influencing departures of sea surface elevation from a level of constant 
B ere Implicitly evident In equation (2.11). 

Initially the current field will cause the density field' of sea water to 
adjust to the apparent rotational force. The density field adjustment can be 
evaluated using hydrographic data through the dynamic height method (Fomin, 1964; 
Neuman, 1968). In addition to geoatrophlc adjustment, surface slope can also 
change In response to gravitational anomalies, l.e. , local changes In the value 
of the gravitational acceleration function. Gravitational anomalies can become 
especially significant for two reasons. Firstly, we are interested In obtaining 
Information about Che current systems based on surface slope measurements only 
end therefore any departures from a flat surface arc interpreted as having been 
caused by geostrophlc adjustment (or winds, etc.) If no other correction has been 
made. Secondly, the Culf Strcom betveen Florida and Cape Hatte^'*? oAnpraily 
follows the continental shelf break countour which Is a region of rather abrupt 
changes In the gravitational acceleration value and subsequently In the gravimetric 
geoldal height. Ac the continental shelf break, there are changes in depth of 
Che order of a hundred meters to five thousand meters over horizontal distances 
of the order of several to several tens of kilometers. 

These two facts, l.e., the existence of Intense current systems along the 
shelf break of western oceanic boundaries, in the Northern hemisphere, coupled 
%rlth the dr.'tmaclc change In gravimetric geoldal height transverse to the western 
boundary make for extremely difficult complexities In the solution of the problem 
of interest to the ocean dynamlclst as It is coupled to the primary problem of 
interest of the gravity minded nvnrlnc geophysicist. Both problems, that of sea 
surface slope associated with winds, currents and hydrography and that of 
determining the crustal and subcrustal structure of the earth can be attacked 
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nose directly and synoptically through the use of satellite imagery versus the 
laborious, time consuming, unsynoptlc (without the co-ordinated use of several 
ocean research vessels) direct sea surface, single ship effort. Unfortunately 
the problems become coupled in actual satellite observations so that initial 
Interpretation and future decoupling of the data are problems which must be 
addressed directly. 

Sea level also varies along continental boundaries as well os transverse 
to the margins. Along the eastern Atlantic const of the U.S., sea surface 
height drops approximately one-half meter between the latitudes 20*N to A0*N. 

It has been demonstrated (Sturges, 1974) that this meridional slope along the 
Inshore lateral boundary of the Culf Stream has the correct sign and magnitude 
to be balanced by the cross-stream gradient of the alongstream current. This 
is an extremely important finding concerning the dynamics of the Gulf Stream 
and furthermore may help to explain why sea level on the Pacific coast stands 
higher than chat on the Atlantic coast. 

Corrections for changes in Che gravitational potential can result in sea 
surface height differences of several meters across the Culf Stream as a function 
of its path. Because of the afore-mentioned complications, initial surface slope 
calculations will have to be based on hydrographic, current and wind data and 
adjusted with geophysical gravity data including Free-air, Bougucr corrections 
and a variety of Isocatlc and hydrostatic compensation adjustments. 

The change of the surface height across the Culf Stream may be of the order 
of one to several meters resulting in a cross-stream slope of at most 10 This 
gradient is too small to be realized by conventional oceanographic methods. 
Consequently, previous to the Introduction of remote sensing, the only sea surface 
work was done with tide gauges and hydrography. The data from tide gauges is 
direct but of itself yields only longshore slope. The hydrography data, though 




difficult to obtain can be used, with appropriate reservation, to indirectly 
compute of sea currents and sea surface slope. This method was generated by 
Helland-Hanscn (1903) on the basis of the Bjerknes circulation theorem (Bjerknes, 
1900) and nukes it possible to indicate surface topography without dire'* 
measurement. The technique was, in its initiation, accepted as truth ’^ut it 
must be pointed out that the so-called dynamic method cannot be used to compute 
* the wind driven current because frictional forces that are not considered in the 
dynamic method, do play a role in the set-up of the pure drift currents. So, 
though the pure drift component of velocity maintains the sea surface slope, it 
cannot be computed from hydrographic methods. 

In a baroclinlc ocean, l.e. , an ocean in which constant pressure surfaces 
and surfaces of constant density are allowed to Intersect, the horizontal 
pressure gradient and current velocity become zero at some depth at and below 
which the baroclinlc and barotropic pressure gradients become equal and opposite 
In magnitude and direction. This level of no motion may extend below the bottom 
so chat the pressure gradients due to Inhomogeneitles in the field of mass are 
not in mutual compensation with the free surface slope pressure gradient. From 
the following equation 






( 2 . 12 ) 


it becomes explicitly obvious that at some depth, at which V ~ 0, the two 
terms on the rlghchand side are equal in magnitude and opposite in direction so 
that mutual compensation occurs. 

The sections to follow explain more fully the theory and actual use of the 
dynamic method. 
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2.2 BASIC cEosTROPtry Arm the PYr^iVMic method 

The opposition of the pressure gradient force and the apparent rotational 
force, the Coriolis force. Is the general dynamic constraint observed by steady 
flow. It is generally acknowledged that steady flow In the ocean has meaning 
only as an average with respect to time, especially ocean currents such as 
the Gulf Stream which observes an extensive ranp.e of frequencies (Pugllster, 

1951, and Werthclm, 1954). Near oceanic boundaries, the goostrophic assumption 
breaks down since sharp horizontal and vertical velocity gradients can require 
lateral and vertical friction to play Important dynamic roles. Still It Is 
possible to assume geostrophy, with informed reservation In a phenomenon such 
as the Gulf Stream. 

Vectorial gravitational acceleration can be expressed In terms of a scaler 
potential function. The difference In the potential function between two 
points lying along the vertical co-ordinate, which Is aligned parallel to the 
direction of gravitational acceleration. Is equivalent to the amount :k 

that It takes to move a unit mass from one point to the other point. . .8 
scaler potential function Is related to the vertical spatial variable by 

dZ« — d«t> » •<> that (2.13) 

% 

levels of constant geopotcntlal correspond to levels of unchanging z. 

In terms of the density Inverse, l.e., the specific volume of tea water, 
the potential function can then be written as 

d4> dp , (2.1A) 

by combining 2.9 and 2.13. Thus the gcopotentlal Is directly determined hy tiie 
pressure and specific volumes functions, which are in multiplicative association. 
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In an environment where all vertical accelerations are Insignificant relative 
to gravitational acceleration. 

In the absence of winds and topographic boundaries, the equations defining 
geostrophlc velocity are 
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+ f -fu = 


ap 

ox 


and 


(2.15) 

( 2 . 16 ) 


which can be written in terms of the potenflal function ^ by Introducing the 
transformation 
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(2.17) 


( 2 . 18 ) 


so that the geostrophlc equations become 
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and 


- 

' OX » 


(2.19) 


+ = -(4^) 


oy 'p 


( 2 . 20 ) 


If equation 2.19 is now Integrated then we obtain 




( 2 . 21 ) 
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which can then be written as 


4>cp) = 4>( p.) - \ o<. dp - \ g dp . 

IN 


vhere o(, is a reference density Inverse of sea water and g is the specific 
volume anomaly. 

Equations 2.19 anJ 2.20 can now be rcvnrlttcn as 


— = - 
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vhere the dynamic height, -ad. Is defined by 


D = f 8 dp , 


if we next consider the geostrophic relations at the ocean surface to be 


- fvio) = - Mssl 


J-Uio') =: - 


a4>(o) 


then equations 2.23 and 2.24 can be rewritten os 
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- -f [ 'U'co) - Arc p) } = _ 


a AD(p) 

ax 


(2.28) 


and 

i } 'U(O) - U(p) \ = - , (2.29) 

J av 

or if the velocity on an isobaric surface, P|^ is known then the velocity at 
any other surface is given by 

■flvCp) -'Vr(p„)) = I ^D(P,)- AD(p) I ^ (2.30) 

The actual. construct of a sea surface picture occurs after the depth of 
the nearly zero gradient current has been determined. Measurements in the region 
of the Gulf Stream vary considerably in depth of actual observation but there 
are sufficient measurements made at depths of two to three thousand meters to 
insure that we've satisfied the criteria of the dynamic method. 

Defant (1941) examined the departures in isobaric surface dynamic depth, 
from a constant value, for a considerable number of station pairs and determined 
that for the Atlantic Ocean, the one thousand to sixteen hundred layer be 
considered the no-motion reference layer. Unfortunately, there is not yet any 
method which applies to all oceans, so one must approach this problem on a case 
by case basis with a belief in relative values. 

Fomin (1964) examined the accuracy of the dynamic method and showed that, 
owing to the accumulation of random errors, computed dynamic heights differed 
from actual isobaric surface dynamic height by errors of the order of twenty 
dynamic millimeters in water depths of a thousand meters off of the Kamchatka 
coast. 
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The accuracy of sea water temperature, salinity and denisty calculations 
^'hould be evaluated, since errors In these measurements give an Initial biasing 
to the results. The extent to which surface topography, as Indicated by dynamic 
methods, departs from true topo is a function not only of the reference surface 
and on the magnitude of horizontal motion at this depth but also, and very 
directly, on the precision of the devices used In actual data collection and 
.iu., . ihis proclr.ion includes the number of vertical and horizontal data 
points, i.e., levels and stations occupied, as well as the more obvious concern 
of instrument scale division. Moreover, the reference surface cannot be located 
at too great a depth because of the Immediate Increase In computational errors 
which results. It should be noted though, that In cases where computational 
errors arc large. It may be of some use to use a shallower reference surface 
in a trade-off of the loss of absolute values versus more exact computation. 
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mometers while salinity Is Implied from chlorine measurements and a chemical 
table (Zubov, 1957). The reading accuracy of deep-water reversing thermometers 
is plus or minus two-hundredth's of a degree celslus while the accuracy of the 
chemical method used to determine salt content Is about plus or minus fourteen 
thousands of a part per mllle. 

It must also be observed that the hydrological profiles across the Gulf 
Stream occasionally extend from abyssal plain ocean depths to continental shelf 
depths, where the depth is less than the reference surface depth. In such a 
case dynamic analyses have been performed (Hclland-Hansen, 193A, Mohn, 1885) by 
replacing the earth's crust of the continental shelf and slope by a fictitious 
water mass and by then extending the ocean density field Into the fictitious 
water mass with artistic license. 
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2.3 Calculations Uslnr, thu Dynnmtc Method 

In order to make a computation of dynamic height and or of geostrophlc 
flow, one must first calculate the anomaly of specific volume, S , which Is 
related to the density of sea water. The density of sea water la determined 
by its temperature (T) , salinity (S) and, to a much lesser degree. Its 
pressure (P) , so that f •• J*(S,T,P). 

Charging the temperature of a parcel of sea water re.'iults In a change In 
the density of sea water by an amount • Similarly is the density change 
due to a change In the salinity within the parcel of sea water. By changing 
both the temperature and the salinity within the parcel of water results In a 
new volumetric density of 




^ P 4T » 


(2.31) 


%rhere the subscript zero denotes the original values of S and T. Note that the 
separate effects of salinity and temperature are not simply superposltlonal in 
nature. There Is also an Interaction term which must be considered. If 

ve now Introduce the parameter pressure, then 


fsTP “ i*S.T.P. + ^ fsp"** 


STP 


(2.32) 


This relation can be reduced by considering the definition of a quantity called 


sigma ~T (0^) , so that 


fjTP = I + 4j>, 1- 


(2.33) 


Using an alternate notation, the specific volume (W) can be written as 
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Syf Sj,p + S 


SP 


>TP SiTp + S^p 


(2.34) 


j’sTP 

and uslns standard notation, a calculation of t» may be done using the reduced, 
or simplified, equation 

0<5Tp =» 0<j5^j^^p-4. AST + S’sp -4- S.fP ^ (2.35) 

where &g.j. is called the "thermostcrlc anomaly." The last three terms contain 
the departure of a from the value ^ p and are known as the anomaly of 
specific volume 


S — Ssp -k- Si 


(2.36) 


>sp OtP ^ST . 

Now variations in 5 indicate the so-called "relative field of mass" within the 
ocean and are used to calculate geostrophlc flow. 

To reiterate, the calculation of <S(z) at two stations must be done before 
the geostrophlc flow can be computed. From the beginning though, it should be 
csphasiccd that derived see hpljhtR and slopes, velocities and transports 

are only arbitrary. The choice of a level-of-no-motion above which the mass 
field and the sea surface are in perfect adjustment, is open to conjecture and 
bc.jmes rather an art based on our present state of knowledge. 

After the values of 5(z) have been computed using hydrographic tables, then 
the procedure for calculating the dynamic height anomalies is as follows: 

First calculate an average S, S • between two depths, say the 

sur a:c, z 0, and another depth, say z h meters, at two different stations. 

< ae can then calculate the dynamic height anomalies (AD) for both stations, where 
in each case AD ■ ^^ 2 * Then the cumulative dynamic height anomaly can be 
assessed by summing the AD values at each station as TAD. Geostrophlc velocities 
and sei surface slopes between stations a and b can be computed. Geostrophlc 
velocities are obtained by the relationship 


^Zi ~ V?, - - 


10 ( ADocC Z, - 2i) - ADfeCZ, - 

f L 


(2.37) 


1 



"T 


1 


■”!' 


1 


1 




15 


where V denotes velocities and L the distance between stations a and b. 

If one wants to Interpret relative flow from Isopycnals, then one could 
consider the formula 


^ Z 


- ^ 
ff 3^ 


( +on y* — t'an 8 ) > where V 


g tan B 


(2.38) 


Equation (2.38) defines a relationship between n^'ostrophlc velocity ami the 

slope of an Isobar (a line of constant density). In the northern hemisphere, 

if the Isobars slope up to the East, then the velocity of the water along the 

Isobar ic surface Is to the North. The steeper the slope, the higher Is the 

velocity. In the case where = o (barotroplc flow) then <'»& and equation 

AZ 

(2.38) may be written In its approximate form 


= _ J 5L tar. ■ (2.39) 

AZ AZ 

Equation (2.39) defines a relationship between the slope of the Isopycnals 
and the velocity shear. For a stably stratified water column, — i->o » in the 
northern hemisphere, the isopycnals slope down to the East when the northern 
component of velocity decreases with depth. Equivalently, we could say that 
the velocity becomes more southerly In character with Increased depth. A 
pictoral representation of this, which represents the Gulf Stream, Is shown In 
Figure 2.1. 

The computer program given in Appendix 1 will now be explained. As 

mentioned above. It Is the convention to compute dynamic dpeth in terms of a 

standard chunk of seawater, in which T • 0*C, S • 35*/,,, plus an anomaly. 

Thus we have the C<_ .,,, term plus the S term. So each parcel of sea 
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water at each depth is composed of two volumes, l.e., each unit volume so at 


.r 



1 


1 


i 

16 


B M 0 meters 


and at 


z •• 100 meters 
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Since the o( p g is a standard and since we are only Interested In differences 

from stations then this program follows the convention of summing only the 
anomalies. 


If we rewrite the anomaly as composed of terms and S 


then 


S = ?T + Ss Sty + S 


tT ^ OttP ^ t>4P ^ o stP • 


(2.40) 


Now In standard oceanographic tables, the first three terms are. combined and the 
last is assumed small, l.e., 


+ s„ , «•“> 

Ss-rp ~ O . (2.42) 

The program works directly with the equation of state and use the relationship 

C(iTP = CX„y„ , , p + S 


Ot JS, o. o * 6. q?2 64- . 

The anomaly In dynamic depth Is then the Integration (2.25) aD • j ^dp . On 
paper the calculation goes something like this: 
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where in the program, the following Is used: 
DIMENSION 

D depth on specific data 

T temperature of data 

S salinity of data 


.1 

I 


I 


Anomaly of dynamic 
cepth ied from 
surfac! down 



1 


18 


-^4 

i 


DD 

TT 

SS 

D1 

T1 

SI 

K1 

P 


100 

101 

102 

1021 

1022 

103 

ICA 


values used by computer In calculating specific volume for 
dynamic calculations* 


values interpolated to standard intervals, i.e., 
0-10-20, etc. 


a counting parameter 
pressure 


DYN 

dynamic meters, ® " j" O^p.T .4 ^ P 

ALPH 

specific volume, p 

NSTA 

station number 

FLA 

latitude 

FLO 

longitude 

ALP 

®^0*C,35V,,,P 

ANOM 

anomaly summation, aD ■ J S dP 

DYA 

dynamic meters, D- j 

DEAN 

specific volume anomaly 

SIGT 

sigma -t 

FORMAT 



data cards 

master card for each station 
printer heading 

H II 

II II 

printout format for depth and velocity 
printer heading 


• r 


1 





■'T 

,1 


1 ■ 


~1 




' I"—'— --T' 


19 


105 

106 
197 
108 


apace 

printer heading 


prlnt-out fonnat for depth, temperature, salinity, sigma -t, dynamic 
meters, anomaly, anomaly sum 


DIN*10 an interpolation Interval in meters — could be changed to a more 
convenient spacing 

11-0 

1. 11-11+1 determines that we start with first station — do the interpolations, 

then look at second station 


1NTERP0L\TT0N 

The Integration of the specific volume anomaly given later in the program 
uses all the observed data, but a value for S and T must be specified at the 
svrfaeo and at standard dcoth intervals. 

This part of the program assigns the shallowest observation to be the 
surface value. Usually this will mean that the first observation at say 0.5 to 
2.0 meters depth will be moved to the surface. Observations between the "surface" 
and the first standard depth are then assigned new index numbers. An interpolated 
T,S value is then found at the first standard depth and so on down the water 
column until observations run out. 

2. Read 100,D(I) ,T(I) ,S(I) 


IP (S(D) 4,4,3 

3. 1 - I + 1 
GO TO 2 

4. MAXI - I - 1 
DD(1) - DI(l) - 0.0 
TT(1) - TI(1) - T(l) 
SS(1) - SI(1) - S(l) 


Reads all the data cards — a blank data card at 
the end gives S - 0, so the computer senses 
that the last card has been read. 


This is the number of data cards. 


This tells the machine that the surface (z - 0) 
values arc equivalent to the first observations. 


T 
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1 
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K - 1 

LI - 2 

DMAS • D(MAX1) 

MINT - DMAX/DIN + 1 
XNINT ■ MINT 
•. IP (•••) 11,12,11 

11. LAST - MINT + 1 
GO TO 13 

12. LAST - NINT 

13. DI (LAST) - D (MAXi) 

TI (LAST) - T (MAXI) 

SI (UST) - S (MAXI) 


counter-index paraneter 
counter-index parameter 
xkaximum observed depth 
number of interpolation Intervals 
real number 


checks for round-off errors 


Assigns to last Interpolated depth the T,S 
values at the deepest observed depth 


Sec the LOxxvwing example, which relates 


to Figure (2.2) for explanation of 


the Do Loops. 


Example : 

Suppose we have the following STD trace: 


S(l) 

S(2) 

S(3) 

S(A) 

S(5) 

S(6) 

S(7) 

S(8) 

S(9) 

S(10) 

S(ll) 


34.1 

D(l) - 1 

34.2 

D(2) - 2 

34.2 

D(3) - 3 

34.4 

D(4) - 4 

34.5 

D(5) - 7 

34.6 

D(6) - 8 

34.8 

D(7) - 9 

35.0 

D(8) - 11 

35.0 

D(9) - 14 

35.0 

D(10) - 15 

0 (blank card) 




•r 
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MAXI ■> 11 *■ 1 <■ 10 data cards 

Did) - 0.0 

Sid) - Sd) - 3A.1 

DMAX ■ D(10) - 15 meters 

NINT - 15/10 +1-2 (integer) 

XNlirr - 2.0 (real) 

DXiAX/DlN + 1.0 - KNIOT - 0.5 (positive) 
LAST - MINT +1-3 
DI(3) - D(MAXl) - 15 
SI (3) - S(MAX1) - 35.0 
DO I - 2,2 
■ XIM - 2-1 - 1 

DINT - XIM* DIN - 1*10 - 10 meters 
MAXI - LI - 10-2 - positive (to 61) 

61. DO 6 J - 2,10 

(D(2) - DINT) - 2-10 - NEGATIVE (to 7) 

7. K-K+1-2 

DD(2) - D(2) 

SS(2) - S(2) (to 6) 

6. J - 3 
continue 
6. J » 8 

D(8) - DINT - 11 - 10 - positive ( to 9) 

9. K - 7 + 1 - 8 

DD(8) - DINT - 10 meters 

SS(8) - ((35.0 - 34.8) (10-9))/ (11-9)) + 34.8 
- ((0.2) (l)/2)) r 34.8 - 34.9 


I 


■r 
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1 


I 


Sl(2) - 3A.9 
Kl(2) - 8 

IF (MAXI - LI) (is 10-8, positive) 
DO 10 J - 8 - 10 
K - 9 

DD(9) - D(8) 

SS(9) - S(8) 

52. KI(l) - 1 
KI(3) - 11 


DYNAMIC CALCULVTIONS 

F(l) " 0.0 eeC sea level pressure - 0 

DEAN(l) ■ 0.0 • let anomaly at surface - 0 

ANCM(l) "0.0 anomaly sum (or vertical Integration) 

ALP(l) ■ ALPHA (0,0,0.0,35.0) is specific volume at the surface 
DYA(l) - 0 

ALPH(l) • ALP1IA(P(1) ,TT(1) ,SS(1)) is total specific volume at the surface 

computed from function ALPHA 

SIGT(l) - (l./ALPH(l) - 1.).10^ is ^ at the surface 

DYN(II.I) - 0.0 dynamic meters 

Do 15 I - 2, LAST 

Do Loops 15 and 16 do a stepwise intcr.ration of the equation of state to obtain 
the anomaly of dynamic depth. 

The specific volume is 

*^P,T,S “ ®^P,0*C,35*/.. + S where % - DEAN(I) 


and dynamic depth 

■> - H.r.s" 


so the anomtly of dynamic depth ap ■ jSdP ■ j dp ^ ^dP -jdp q»q 


• r 
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or 

ANOM(l) - DYNdl.I) - DYA(I) 

2.4 ADVANCED GKOST ROPliY AN D GULF STRL'.AM Tli roUY 

It Is convenient here to Introduce the concepts of baroclinlc and barotropic 
velocities since the former indicates relative geostrophlc currents and the 
field of mass while the latter in related to deep water velocities n.iJ sea 
surface slope. To differentiate fully between the two tcrraonological concepts 
as they apply to a western boundary current such as the Gulf Stream, it is of 
benefit to develop some theory explaining the use and meaning of the concepts. 

To more fully appreciate the nature of the Gulf Stre n it would also be of 
worth to consider not just the dynamics of rotation balancing horizontal pressure 
gradients but also the effects of winds and bottom topogranhy (see section on 
Gulf Stream models) since presumably (sec the section on Results) winds and 
topography do modify the intensity, shape and position of the Gulf ;.ieam. 

If we consider the steady system of equations 

= - -|^ + 

y+u. =- 

and 

- - er «•«) 

trhere Vxz Reynold's stress terms and we've a right-handed cartesian 

co-ordinate system with x cast and y north. 


. r 
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and 


Next we define mass transport as: 

U = f fu.dz 

V = P f U d Z where b 


(2.47) 


(2.48) 


Indicatt'O Lho bottom depth. Then evaluatlu;; the stre;;:. difference 1 rom tc.p.j^, 
to bottom, Zj^, we integrate the central region equations over depth, so we 


obtal-. 


-^v --J 


bP 

ax 


dz + 


(2.49) 


and 


iU 


= - r> 


yz I- 

• '^b 


(2.50) 


where a porti'^n of the velocity components. Individually, balances the pressure 
gradient (gcostrophic part) and a portion of the velocity balances the stress 
(Ekman part). 

Next ve can break the geostrophic flow into both barotropic and baroclinlc 
parts but we must first reformulate the pressure gradient. Consider the 
continuity condition to be 


a fUL 

ax 


a fv . af ^ 




az 


- o 


(2.51) 


Then integrate equation (2. SI) in z to get 


dz = O. (2.52) 


i 
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Nov If VC assume that the volume of water passing through unit width In a 
vertical cross section in unit time is constant with the same condition holding 
for the other two directions, i.c., that 




ATL - 


and 


dX 




]’(Zb)'Ut Zb> 


azb 

ax 




=. o 


o 


(2.53) 




then ve get: 


au ^ aV 


o 


ax ay 

80 there is no mass flux through the bottom or the surface and continuity is 
satisfied by the horizontal mass flux. 

Now, for a re-formation of the pressure gradient, we must investigate 
the term 

n 


(2.55) 


- i ; 


ax 


dz 


(2.56) 


which is equal to 

ap 


-1 


dz = - 

, ax ax 

Zb 


) 


P Ail _ p 


ax 


azb 

*» ax 


(2.57) 


As Is evident, part of the pressure, P, is independent of the density distribution 
and part is not so we must get a reference pressure and get an anomaly from 
that reference where this anomaly does vary with position. Define 


P = 6 j f ^ 


(2.58) 


which is the hydrostatic condition, and 


1 




Ep = j paz , (2.2 

which is the potcatial energy of a column of water of unit horizontal area. 

It thus becomes obvious that we get potential energy anomalies where we have 
pressure .anomalies. Combining (2.58) and (2.59) then yields 


psap , (2‘C 

^ Potm ^ 

where o(. Is a function of pressure (P) but not a function of temperature (T) 
or salinity (S). Note that Is taken at T ■ 0*C, S ” 35*/»* or 0 ( 35 ^ 0 ,? 
and S is the specific volume anomaly. With these we can go back and rewrite 


equation 2.57 as 


,Pk , 


_ I ' i£di - ^ ( '4- <*p 
J ax j c 




pe _ >Zk 

“5“ ^ P “ dx 


where this Is 
a function of 
bottom pressure 
only 


where the integral 
can be called the 
potential energy 
anomaly, i>. 


where “Of so we obtain 


X dx e>x , . 


- P, 


» ^ 2b 

ax 


i.e. , a<, - ®35^o,P 

Is a function 
of bottom 
pressure 

Now, ^ Is the only term dependent on the T,S distribution directly although 

ax 

P. (bottom pressure) does depend on water density above the bottom as well, 
b 

We can now combine terms to get 

- ^ _ p / 3Pb , aZb V (2.63: 

J av ax ' 9 ax ' • 


represents the 
barocllnlc part 


represents the 
barotroplc part 


I 
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where both parts of the barotroplc term arc related to the bottom slope. 

«>Pb 

is the bottom slope directly and is equivalent to slope since it Is the 

change of height of the water column above the button. The terms are essentially 

jf opposite sign since Pj^ increases down and z increases up. These arr both big 

terms but the difference between them may be small. If we now asF'imo tha": 

z C3i z. , l.e., a reference level near the bottom, then we’ll be able to get rid 
b 

of a largo, part of the.- term. Recall that, 


Pb = Pz + 3 C2.64) 

where the first term is written as a reference pressure term near the bottom. 
Thus, 



( 2 . 65 ) 


wnuLu 


XXU14.U Ul. XUCC;^L«&Ci.OiA ^ Xd xttVAC|^wttwv..it u wx a. ^ aw X 


dpb 


ex 


= 

ex 


2= Zb 

If we now multiply through by 
we obtain 


z «. Zb 
e.».lPo) 


- ir • 


(2.66) 


and substitute the result into 


]' 


ex ex 


p.( 


0 *.(Pb) 3P2 


3 


Z ■ Zb 




( 2 . 67 ) 


where if we have a bottom where T and S arc essentially uniform and use this 
as a reference, i.e., take at the S,T of the bottom of the ocean, then the 

last term on the right cancels out since o<,lP,) = . Unfortunately, this 

Jb 

decoupling breaks down where there are significant intrusions of deep water 
near the surface. Finally we obtain 


I 


J 
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Pb a Pi 

9fb 


( 2 . 68 ) 


barocllnlc barotroplc 

part part 

JP 

where — - is depth independent. If integrated, ’ this indicates a uniform 
velocity from top to bottom with a ir.o^.nltude equal to that at the bottom as shown 
in Figure 2.3. We can new write the Integrated equations as: 


the Ekman transport relations, 


“ j- — 'C’ ! 


2-1 


(2.69) 


4 Ue = '?:y = 


2 » 


the Barocllnlc Geostrophlc relations. 


- ^ 


(2.70) 


(2.71) 


i-Uq = - 


(2.72) 


and the barotroplc geostrophlc transport as 


- fVB = 


P b 3 ^ 
9^6 


(2.73) 


4 Ub = - — 

6f, ay 


(2.7A) 


where the subscripts B,g and E denote barocllnlc, barotroplc, and Ekman 
transport, respectively. 
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Summed up, we have a velocity satisfying the total relation 

U = Ut Ug -*• Ub. (2,75) 

The 2-D continuity equation is not satisfied by any of these velocities 
separately except in special cases. Rather, In general, the continuity 
condition is 


a \J ^ av 
ax ay 


o 


(2.76) 


In the case of Barocllnlc Geostrophic Transport (BCLT) , there is a divergence 
because of the Coriolis parameter's dependence on latitude. In the case of 
Barotropic Geostrophic Transport (BTT) , there can be a divergence due to changes 
In depth. So we can get a number of variations depending on how we drive this 
with vln'i We can drive it such that the wind driven part satisfies 

continuity alone or such that we need the geostrophic part, i.e., if Vh TJ e == O 
nothing happens but if then we get the balance from other components 

(geostrophic) . 

The break up of motion into such parts is not quite complete but it includes 
vhat we use in ocean circulation studies. 

There are exceptional cases where a bottom Ekman is important. This 
occurs where a bottom stress causes a significant transport in the bottom layer 
and occurs in cases where Barotropic flow is dominant over Barocllnlc. 

We can have the motion composed purely of one component when the divergence 
of that component is zero. Otherwise, we need enough of the other components 
to satisfy the continuity balance for the total motion. 

Please note that in the open ocean, the only way to get an Ekman transport 
1 b to have a non-zero curl of the wind stress. 


® Note that f will be assumed a function of latitude as f " f 
where ^lud 0 aie constants ( p* 10“^^). 


local 


+ By. 


r 
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Thus, by cross-dlf f crcntlatinB the component equations, we can form 
vortlclty relations which contain terms that are the horizontal divergence of 
the particular component of motion. So, for the Elcman Transport, we have 


f ( ? . ^ L + , + a y, ~ 

^ ^ ^ 3 K ~ dy 

> J t.. . t t— . I. , - . . ■ — t 

<a) tb> I c ) 

where (a) is the rate of change of vortlclty due to stretching, 

(b) is the advcctlve change of vortlclty, and 

(c) Is the production of vortlclty due to torque (or curl) 
of the wind stress applied to the sea surface. 

This equation can also be written as: 


llh 4. 

dy 


d ( Tyz 3 t ^ 

^ ~~Y~ ) ay ^ ^ 


^ horizontal divergence 


■ ^curl of 


/Wind stressv 
V f ) t 


( 2 . 77 ) 


( 2 . 78 ) 


80 we can have the Ekman velocity in pure form only when the curl of 
^wlnd stress ^ ^ q 

By cross-differentiating the ECLT we obtain 


i ( ^^9 

^ ' 3 H 


?.y^ ■■ \ 

ay > 


f Vg 


o 


( 2 . 79 ) 


which says that the horizontal divergence is equal to zero only when there is: 

( 1 ) no meridional velocity component, and 

( 2 ) no dependence of the rotational parameter on latitude. 

The curl of the BTT yields 



V.. , I « _! t 


ta> 


tb) 





i 
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whore (a) and (b) must be considered non-zero terns since we may not have a 
zero chance in depth and since bottom pressure can change significantly with 
position. In (a), is proportional to and in (b) , is proportional 

to so the right-hand side ((a) + (b)) represents an advective change. 

Thi . Barotropir. Geostrophlc balance can also be written as: 

^ + ^ _f ( y, - Vg , 

^ ^ \ m I ■ I.I. . -I ..4 

la) cb) CC) 

where (a) is stretching due to advectlon, 

(b) is stretching due to local change of earth's rotation, and 

(c) is stretching due to change of total depth, 

so the horizontal divergence is zero when either there is no meridional component 
and no change of f with latitude and no change in depth. 

Ekman transport takes place within a certain distance from the boundary which 
is essentially the Ekman depth. This depth is dependent on the eddy viscosity 
and the coriolis parameter. 

Tlie boundary generates shear waves which propagate into the fluid. This is 
the same kind of boundary as with an oscillating plate — but the boundary is 
rotating versus oscillating. So the velocity decays with depth. The problem 
is made time independent by eliminating the time dependence on the rotating 
plane where the time dependence really is. Note that we can obtain the 
meridional component of velocity directly from the wind stress curl. This is 
the vertical vortlcity balance with no stretching. 

If there is zero wind stress curl then we have only one velocity component 
which is in the x-directlon and will be a function of y only. Since if the 
wind stress curl is zero then Vg ■ 0, which indicates that U • U(y) only so 


r 


I 
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■=o» Thus VC get a zonal flow varying with latitude for zero wind stress 


curl in the pure Ekman case. 

For the above we hsivc no boundary effects. When we have boundaries, then 

this type of notion (flow) doesn't satisfy the boundary conditions so we may 

need another type of flow which is also dlvergcnceless in the interior regions 

to satisfy the boundary conditions. Such a flow is called coastal ujiwelling. 

We can have pure baroclinic only in cases of pure zonal flow, l.e., flow 

only along contours of constant coriolis parameter f. If V >0 then U f 

S S 

\> but U “ U / V , But when we introduce boundaries, then this is no 
g(x)* g g(y) only. 

longer possible. Pure baroclinic flow is valid in the atmosphere, the Interior 
of the ocean and in zonal canals (such as the Antarctic Ocean) . 

Then, for pure geostrophic barotropic transport. 


-i— ( i- Ue 'I - ^ f J- 
' h “ ' ay ' h 


Vo 1 = o 

- / 


(2.81) 


and the flow must be along lines (contours) of (f/h), otherwise the flow is an 
arbitrary function of (f/h), so if the bottom is level then the flow is zonal 
(along constant f) and if f doesn't change much, then the flow follows the 
bottom contours. Since (f/h) is constant along a streamline, for this case, 
then the flow goes towards the equator in crossing a ridge and towards the North 
or South Poles in crossing a trough in the ocean. 

We can also have combinations of ET, BCLT and BTT such as: Stress-driven 

Barocllnc Transport. 

This is the type of problem that is considered for most ocean circulation 
studies. We h.avc both the ET and the BCLT systems of equations which are 
combined to yield the Sverdrup Equation: 


p t V, . V, ) = ^ 


• XZ 


. r 


(2.82) 
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So, for Che Sverdrup Equation, wc require Che condition Chet 

aUe ^ aVc _ 

ax ay 



vhlch Is the continuity equation for the total flow. 

a "Cvz 


P ( Vf . - ) 


ax 




(2.83) 


(2.84) 


which is known as the Sverdrup Equation. 

There is also Che possibility of an ET and BIT type of flow represented 
by the following balance 


PV. - ( 2^ 


2^) = - hCUe|j(i)- . (2.85) 

k_ I 

the change in vorticity from the 
barotropic position due to scream-> 
lines crossing contours of (f/h). 


There is a limit on doing this kind of flow because in most oceans bottom 
friction may be small-since most of the current is normally picked up in the 
baroclinic portion. But, neglecting baroclinic flow does require the inclusion 
of a bottom stress. 

In applying these concepts to actual Gulf Stream models, it is proper to 
consider first, Sverdrup's Theory of Wind Generated Transport. The equations 
of motion governing the system are 


and 




( 2 . 86 ) 


(2.87) 


which ore integrated from depth z ■ -h to the surface, where 7 . " ^ 0. It 


is next assumed that the stress at z ■ -h is zero, in accordance with the 





( 2 . 88 ) 


( 2 . 89 ) 


whorc' the* subscii.pt: s linpliLS the f car.ur f-'ice value of the wliul stress. Then 
taking the curl of the equ.itiuns of motion and adding the resultant rcl.itionu 
with the continuity condition, the following results 





( 2 . 90 ) 


■ which says that the North-South transport is completely determined by the curl 
of the wind stress. If this relationship is now substituted into the continuity 
equation and an IntcgLULxuu wlLh Luapex.1 to x la pc- rforiucd^ wHcn 

^ 5 -fy ( )z 

where g(y) Is a constant of integration, now g(y) is an arbitrary function of 
y which must be evaluated. The equations are of first order so we can satisfy 
only one boundary condition but by the nature of observations, we know that in 
the western part of the Atlantic Ocean, there is a narrow band of intense 
current occurring, so this solution will probably not apply to the western boundary 
region, or at least the distance c .away from the boundary which is at x - 0. 

Though this early model of wind generated transport seems inappropriate to 
a Gulf Stream model, it did establish a theoretical basis from which to proceed. 

So we've solved a circulation problem with U being zero at one coast 
(Eastern) but we were unable to satisfy a condition on tlie other coast (Western). 
The dyn.amics of the problem were not of a high enough order to permit a clcscd 


!■ 
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basin solution. 

What we really have here is a divergence of the F.kmnn wind drift, produced 
by the convergence of the gcostrophic flow (or vice versa convergence by 
divergence). Thus the picture is of a coupled flow comprised of the Ekman 
layer responding to the wind and the underlying gcostrophic layer which adjusts 
Itself. 

We'll next look at a "closed" circulation In a bounded basin— son.othlnR not 
possible in the bverdrup approxicutlons. Stommcl (1948) introduced a frictional 
Cera or rather, a mechanism for dissipation. The friction is introduced in 
Che form of a drag on the bottom of the ocean — which may or may not be a 
realistic dissipation source — to show that closed basin circulation with 
western IntersificaClon arc possible. Consider a rectangular basin, Oixi L, otviw, 
and a wind stress given by Xty) = - F We now have 


” h '”n 


(2.92) 


f -I t I A* = (Ax^)j'' =Tv>-T:y'> 

I L. 


(2.93) 


If we next choose 


rcv) 

f 


= - F ) 


f 


= RU 


U 

-h 
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R. V 


AX m 



AX d2 


«nd R is a constant measuring the bottom drae. Stommcl then integrated from 
* • -h to 7 ■ using, by definition; 

h aI 

5 *^ ATdz = V hv 

-h 


and noting that 
'h 


4) ^ f 4 ^ ^ ° 4 ^ 


c 

in 


= - 9t^ - - s^-f? 


And 


- _££. dz — L ae 

-h^ f Ih ^ 


lldz 








so the main equations become: 


i Vi A) 

= - 9h 

4 [ - F tosC"-Vb>] - 1^^ 

jy\u 

=s - gh 

ia - V 
»y 




dx 


- •^r 


(2.94) 

(2.95) 

(2.96) 


Taking the curl of the equations of motion and summing yields in terms of a 
stream functlon^^ such that 



V “ 




'U 




results ir the following equation 

hB a'l'' uF , 3*'^' . 3*'^' 


o 


to be solved. 

If It Is now assumad that the y-depcndcnce for ^ is of the form sin 
and If we assume that ''V' can be written as t^^/b) 
then substituting this Into the equation yields 




dCr / Tt 


d X 


- , 


wnulu 


Ot 


h b 
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Ond 


V' ■= 


TIF 

bR, 


Now, the general solution to the homogeneous equation 2.97, with 


A - - 


o( 


ot' 


TC* '/i 


r)" 


and 


Is 


B- - (-£^4 

Z ^4 ^ 

Cr ~ Cl C2 e®’' 

where Cj^ and C 2 are the constants of iTitcgratlon , so 


'M/'u.'i) = •+. - I ) ^ ) 


icy 


Ca - 





(2.97) 




(2.98) 


(2.99) 


where 
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For the equation, the b^^undary conditions arc that 

“0 at y ■ 0, b 

and X * 0, L 


so that 


''4’^<0,y) - ^ (h.y) - A.1^' (x,0) - At'^(x,b) - 0. 


The condition on y “ 0, b Ir 'it immediately since the result is a unction of 
sin . The condition on x ■ 0, L yields that 

p 


Ca = 


l - & 

e®*- 


and Cg = 


e*" - . 
e*'- - e®'- 


SO the Stonunel solution is that 


I e'-T"' -a-' ■*-?'>''• I- ) - 1 

I ( , . e'- ^ ‘ ) e ^ ^ 




and thus 


-r--^ I ae*'* Cse®'- , ] C04 ^ 


V = ^ -^.lA c. 6C, e®" ] i.'o 


In Interpreting this, Stommel chose (rather arbitrarily, to say the 


least) , 


L ■ 10 cm 


b - 2 (1C°) cm 
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h - 2(10^) cm 

R - 0.02 

2 

F ■ 1 dync/cm 

Nov If the ocean Is not rotating, l.e.. 




P = 


then 


ana 


C4 = e 


- ’^yb 


-M-' - =■'" ( ^ ) ( e 


' b 

a solution depleted In Figure 2.4 
The sea surface Is given by 


-S’ yre. 

-i- e o 
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( 2 . 102 ) 


and depleted In Figure 2.5. Now when the Coriolis parameter is a constant 
(0 .25 (10 ^)) then the streamline diagram remains the same as Che ^ *• 0 

case but the diagram changes as shown in Figure 2.6. Next, when the Coriolis 
parameter Is a linear function of latitude, i.e., f “ fo ^ y then Figure 

2.7 results for and Figure 2.8 is the sea surface topography depiction. 

The obvious feature of the f constant case is the crowding of strciirallnes 
on the western side of the basin. The velocity along the Southern border is 
^20 cm/sec and along the Western border is 240 cm/sec. The band width of 
the Intensified current is 100 km. The Implication is then, because f changes 
as a function of latitude, then there is a concentration of streamlines along 
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the western boundaries of oceans ultlnR In such phenomena as the Gulf Stream 
and tha Kuroshio. 

We'll now consider the Munk (1950) frictional theory of the ocean with a 
model structure similar to that of Stommcl's. Hunk's motivation to do this 
problem was that Stoiumcl made the assumption of a non-zero drag at the bottom, 
which might be ratln^r unreallr.f ic , so Munk solved a problem using the assumptions 

that at some depth z = -z , the currents vanish and the stress Is ascumevl 

nor;o 

negligible and that lateral friction is not negligible. The equations of motion 

are now 
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Munk then used the concept of vertically integrated mass transport which 
was introduced by Sverdrup. After Integrating in z, cross-differentiating in 


X and y and combining the equations in • where V 
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then the equation to bo solved is 
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for singular boundary region 


for the interior region 


This equation is now applied to the rectangular ocean. 

The boundary constraints arc that both 'Vp* and its derivatives normal ^o 
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the boundaries oust varnish. Note that we've assuocd only an cast-west wind 
system 60 

Hunk's rational was to look for a boundary layer solution of 2.106. 

The motivation to do this comes from the problem's physics which Implies that 
friction may play a small role In the interior of the ocean so In the interior 
Kx may be small. 

From Figure 2.9 , Txs can be represented by 


Txs = -T: cob 


(2.107) 


A change in variables such that 

y' ^JL 
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is now introduced. 
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Then if we assume that is of the form Smay^as a V* variation, then 
/ 

we can write the x -dependent oart of 2.106 as 
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and furthermore, if we call ^ ^ the Interior solution then ^ C . 

Near the . undary, frictional terms are Important and since the boundary layer 


is small we can stretch the x co-ordinate in this region by introducing the 


variable ^ defined as x'=x"’^ where m is a constant determined from 2.108. 
Since we require the balance In the boundary region to retain the — - 


term. 


then we would like this term to be of order unity in the boundary region. We 
now match this term with all others in the equation which arc of order unity. 

am,* , . . 

If this is done then we find that and balance each other fairly 
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veil. Equation (2.108) then becomes 
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and the balance in 2.109 Implies that 
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and to an order of , this equation can be approximated by 
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for the boundary layer. 

The final solution of 2.111 Is 
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and if we apply the boundary conditions that Aj/^o on all boundaries and that 
bAb' 




■••oat x= 0, a then the complete solution is 
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So using the mean annual zonal winds only, the solution shows that the 
integrated oceanic wind-driven circulation is divided into closed circulatory 
systems called "gyres.” The gyres arc bounded at latitudes where the curl'^T = o 
and are centered at latitudes where curlT is an extremum. 

So for the zonal wind case given here, the dividing lines between the gyres 
are given by 
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and the maxima or the axes of the gyres arc at 
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The form of the stvcamllncs is as shown In Figure 2.10. If we define 
K ■ ( then the region between x •• 0 ami x “ A/K corresponds to the 

Gulf Stream. Now If the solution that we have Is written as 
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then F(Xj^) Is the total northward transport between x - 0 and x - x, and 
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Is the total northward transport per unit width between x ■ 0 and x ■ x. Now 
dF 

when F and -r~ arc evaluated then wo find that the ocean falls Into the three 
dx 

parts (as a function of x) which were used In the solution development, so at 
the western boundary then x << r and r » 1, and for any fixed latitude, 
varies with F(x) only, so 
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This I'cprcscnts an undcrdampotl oscillation with wavelength 
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Note also that the form of results in a counter-current outside of the 

7 2 

western boundary region and for 5(10 )cm /sec L 200-250 km. 

Recent data suggests that the Gulf Stream is 50-60 km wide, so 10 

is a more reasonable characteristic lateral frictional transfer coefficient 
magnitude. 

Away from cither the east-west boundary, x is comparable to r and 
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which give a broad constant drift 


which corresponds to the Sverdrup solution where 
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80 the width is small, the current is weak and in the real world a local meridional 
wind may obscure the current. 
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Hunk found Chat the countercurrent cast of the Gulf Stream was In good 
agreeaent with theory but there Is also a counterflow observed also between the 
coast and the strong meridional current which Is not predicted by Munk's viscous 
boundary layer theory. 

Although the frictional baroclinlc models are quite good there are problems 
such as: it 1s assumed that sufficient potential energy is available to allow 

the circulation to exist as a baroclinlc mode, but it i. ly not be Lruo that for 
every distribution of wind stress, the density field will respond with a 
distribution of potential energy that is determined by Che wind stress. 

Since the non-linear terms are neglected, friction acts the same up and 
down stream along the western boundary, and the question remains "why does the 
western boundary current separate from the coast." 
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2.5 NU>tKRICAI. MOni-I.S 

Though the analytical Munk (1950) node! successfully predicts, to varying 
Measure, the gross stationary features of the Gulf Stream, it Is evident from 
ocean station data as well as present VilRR (Very High Resolution Radiometer) 
data that such a model is a simple approximation at best. There have been a 
'• multitude of models either in extension of the Munk work or with non-linear 
approaches to the problem, such as those of FofonoCf (195A), Morgan (1036), 

Charney (1955) and Stommel (1961), but no one has yet constructed a theoretical 
model in which one can have more Chan reserved confidence. 

In all of the above models, non-llnearltlcs are not considered contemporaneously 
with lateral and or vertical friction. Furthermore, temporal variations and baro- 
clinicicy are neglected entirely. It thus becomes implicit, from the lack of 
good reproduction of nature by theory that the stationary barotropic assumption Is 
not a very complete form of truth in approach to an investigation of the Gulf Stream 
nature. 

It Wes, in fact, not until Lineykln's (1957) work appeared that the "diffusion 
of density" equations appeared in any of the mathematics systems approaches to 
the problem of solving the complete system of coupled equations (Nos. 2.1 through 

2.7). 

Sarkisyan (1969) indicates the more obvious shortcomings of oarotropic 
oceanic circulation models. He demonstrates that failure to incluac baroclinlcl ty 
and bottom topography arc perhaps the most serious omissions. It is clear from 
his discussion that the density field, more than any other variable, defines the 
current field, so flow-field calculations are only as accurate as our knowledge 
of the field of mass. Consequently, Bryan and Cox (1967), Bryan and Cox (1968), 

Bryan (1969), Sarkisyan and Keondzhyan (1972), and Kcondzhyan (1973) took the 
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numerical methoJa approach to calculating sea level, mass transport and flow 
fields using numerical, versus purely analytical, techniques. The latter two 
works furthered the numerical predictions with a diagnostic approach to the 
problems. 

The present state of the numerical modelling art is not sufficient in terms 
of computer storage space, time limits, mesh grid and finances, to solve 
equations 2.1 ll.rough 2.7. Presumably with the advent of larger and fanter 
computers, we'll be able to more adequately predict tlie velocity and deusily 
fields but presently the straightforward numerical approach has glaring short- 
comings . 

The diagnostic approach to numerical modeling prediction has been taken by 
several investigators based on the realization that there is a considerable 
amount of ocean station data which has been collected over the last sixty years 
in the near vicinity of the Gulf Stream. Given this data, Kcondzhyan has 
reconstructed the density fields as a function of time and then directly computed 
the velocity and pressure fields which arc contemporary with the STD field 
observables. With such an approach, one is not burdened with non-linearities 
in the equations of salt and heat and in the non-linear coupling of S and T 
with the horizontal equations of motion by way of the quasl-hydrostatic assumption. 
Keondzhyan notes that the results can be assumed to be a multilayered model, , 
which in the limit of zero layer thickness becomes a representation of a continuously 
stratified fluid. This approach is further supported also since no layer — inter- 
facial conditions need be imposed other than the continuity condition. 

In the actual model, Kcondzhyan imposes a wind stress and supposes the 
conditions of no flow normal or tangential to the lateral or vertical boundaries. 

With a basin defined laterally as and on the bottom by M(O.x) , where 

® and A are latitude and longitude respectively, the author then specifics the 
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density field Q,S,Z) nnd nuncrlcally solves a modified version of the syatc^n 
of hydrodynamic equations 2.1 through 2.3 with a horizontal step spacing of 1.25". 
The calculation results yield velocity profiles which are in reason.iblo agreemont 
with present knowledge of the multlyc.ar averages for the summer season for the 
direction and strength of circulation in the North Atlantic, The cyclonic nature 
of the lowest order circulation i.s evident with the East Greenland and l.aborador 
currents and a strong Gulf Stream ..howlng up quite clearly. Flgtire 2.11 depict', 
the field surface elov.-itiou as covnputed froai the predicted field o£ pressure 
anomaly of Kcondzhyan. The predicted change in sea surface elevation across the 
Gulf Stream is then used In a one by one comparison of surface elevations computed 
by the dynamic height methods from actual station data, as discussed earlier, 
along identical transects of the Gulf Stream. These comparative results are 
shown in all figures entitled Sea Level Elevations in this report. 

A startling fact arises from these results, as evidenced in the figures, the 
Kcondzhyan (K) numerical model, though only useful in a time-averaged sense, 
shows remarkable resemblcnce to the dynamic height calculations generated by 
Huang and Pietrafesa (H-P) in this report. The curves are never identical and 
the departures of the 11-P curves from the K curve are not well correlated to 
the contemporaneous wind magnitudes and directions. It must be remembered here 
that the K curve includes a seasonally averaged surface wind so that departures 
of the mechanical forcing function from the mean should be the rule rather than 
the exception. 

The results of the Kcondzhyan work depict the following features of the 
North Atlantic ocean: the well-known counter-clockwise Northern cell; the East 

Greenland and a substantial Iwiborador current; the clockwise rotational circulation 
at mid-latitudes; an Intense western boundary current, the Gulf Stream, which 
separates at approximately 50* west longitude and 40* north latitude into a 
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northeastward flowing branchy a continuation of the Culf Stream and a part which 
turns, In a clockwise fashion to the south, to become a current counter to and 
beyond the eastern lateral boundary of the Culf Stream. 

The numerical models also substantiate the barocllnlc ocean work of 
Llneykin (1955a, 1935b, 1955c, 1956, 1957a, 1957b, 1957c) In contrast to the 
classical theories based on the Sverdrup relation. The results of Kctmdzhyan 
also show, rather conclusively, that the meridional and zonal transiiorir are 
not independent of one or the other, as Is suggested by the Sverdrup theory, 
and extensions thereof, but depend. In a coupled way, on wind stress, bottom 
topography and the barocllnlcity of the ocean regime. This diagnostic approach 
offers a powerful tool for future work In ocean dynamics. 
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Figure 2.2 Salinity versus depth for cxiimplc given in Section 2.3. 
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Figure 2.6 Contours of sea surface height 
in cm. In a uniformly rotating 
ocean (Stommel, 1948, Fig. 4). 



Figure 2.7 Streamlines for case in which 
Coriolis pnram'tcr varies 
linearly as a function of 
latitude (Stommel, 1948, Fig. 5). 
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Figure 2.8 Contours of sea surface height 
in cm. in an ocean where the 
Coriolis parameter is a linear 
function of latitude (Stommel, 
19A8, Fig. 6). 
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Figure 2.9 Rectangular ocean and 

wind stress, .as applied 
lat i t udina 1 1 y over the 
ocean, in Hunk (1950) 
model. 
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Figure 2.10 Mass transport strcamfunction 
in units of 10^2 pram sec“^ 
for the Munk (1950) model of 
the Gulf Stream (Munk, 1150, 
Fig. 5). 
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Figure 2.11 Theoretical Model of Gulf Stream Surface 
Elevation by Kcondzhyan (1973). 





3. SURFACE PROFILES MSllD Ori HYDIUX^nALHIC DATA 


As discussed in the previous chapter, the surface elevation of the ocean 
can be related to water movetnont direct!;,- throuj'h rlic geostrojihic equations. 
Accepting the approximations of the inod< I. we hove calculated the ocean surface 
elevation profiles based on hydrographic data tabulated by National Oceanographic 
Data Center (NODC) , NOAA alonf* the Gulf Stream path off of Che eastern const 
of the United Statco, Tha u .tr. .va t • ' ’>/ , ->h?c covering the 

period from 1912 to 1972. The location eta. . on i... i-i.r -iC plotted on a 

nap with time of the survey indicated so that the ship track of each individual 
cruise can be identified. All of the stations taken on a particular cruise 
are then fed into the computer as described In Chapter 2 so that dynamic 
heights can be calculated. 

From the station data, one finds that there are many cruises designed to 
study continental shell water only. On such cruises, ship tracks stop as soon 
as the edge of the Gulf Stream is reached. Since our interest is, however, 
in the location as well as the magnltudi> of the surface elevation deviation 
caused by the Gulf Stream, any ship track that duos not completely traverse 
the Gulf Stream is discarded. Thus the data from a number of very closely 
spaced stations collected mainly by Duke University (Stefannson, Atkinson, 

1967) is unfortunately, of no use in this study. Also, since one of the basic 
assumptions of gcostrophlc modeling is the existence of a layer of no motion, 
models based solely on shallow water data will become increasingly Inaccurate 
when the depth becomes too shallow. As a result, we exclude nil of the stations 
with a total depth loss than 30 meters. As for the layer of no motion, an 
arbitrary depth of 1000® is used throughout the present study. When the 
total water depth is less than lOOO"*, the bottom is used as the layer of no 
'motion. This may not be the best choice, but it is adopted as a compromise 
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BO thac the vast amount of data «.ollcctcd by different ships at different 
times and locations can have a uniform reference. 

All of the results are presented in graphs on a one for one graph to 
cruise basis. Sometimes when a single cruise covers two tracks, the profiles 
of both tracks are plotted together to show the difference of surface elevation 
over the short time period between the two transects. Wind data arc presented 
uoi^cther with the profile:: whenever available. The wind :;pcod Is measured In 
knots and the direction in degrees with zero pointing to the north. Together 
with most of the profiles, Keondzhyan (1973), mean Gulf Stream model results 
are plotted as a common reference. The results are grouped by geographic 
locations and are presented as follows. 

Figure 3.1 Indicates the ship tracks of the survey conducted in 1932. 
Figures 3.2, 3.3, 3.4, 3.5 are the earliest complete profiles available from 
the NODC data. The ship tracks followed the line between Cape Charles and 
Bermuda. The profiles show considerable fluctuations In Gulf Stream strength, 
but the location of the Gulf Stream is rather stable at 200 nautical miles off 
of the coast. Unfortunately no wind data are available for these cruises. 

Figure 3.6 indicates another eight ship tracks In the same general area 
off the coast of Virginia covering the years 1953 to 1966. Again, the agreement 
between the profiles based on the hydrographic data and the mean Gulf Stream 
model of Keondzhyan Is acceptable. In figures 3.8, 3.11, 3.13 the ship tracks 
started and ended In the middle of the ocean. An arbitrary station in each 
cruise Is picked as the reference station. The distance from land marks to 
such reference stations Is indicated on the graph. 

Figure 3.15 indicates ship tracks of all the stations south of Cape 
Hattcras covering the years from 1955 to 197C. FigureR 3.16 to 3.30 show the 
profiles, the available wind conditions and the profiles from Keondzhyan's 


mean model. 
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Of special interest are the stations indicated in Flp^urc 3.31. These are 
rather closely spaced stations both in the time and space domains. Figures 
3.32 to 3.36 show the ten profiles surveyed in eleven days in June, 1968. 
According to the accepted theory, the Gulf Stream path will not change much 
in such a short time span as this (see, foe example, Stommcl, 1966). However, 
the hydrographic data showed that dramatic changes could occur over a very 
short time. Figure 3.32 presents two different profiles obtained over the same 
transect within a three day period. The remarhable thing is that surface 
elevation increased by half an order of magnitude from one curve to the other 
leading to inevitable questions concerning the accuracy of shipbound surveys 
spanning months during a single cruise. Similar features of short period 
changes have also been observed by recent remote sensing techniques. Such 
phenomena undoubtedly deserve further scrutiny. 

Another group of closely spaced station cruises are presented in Figures 
3.37 to 3.45. There are nine cruises transecting 17 tracks surveyed from 
August to December of 1965. Though the profiles in general agree very well 
during any single cruise, the differences between cruises are still substantial. 
This again indicates the mobility and transient nature of the Gulf Stream even 
south of Cape Hattcras. 

In general, there seems to be no direct correlation between the local wind 
condition and the Gulf Stream strength. This is supported by Stommcl 's 
(Stommcl, 1966) energy argument. This leaves rapid changes in Gulf Stream 
Strength still not accountable. Another interesting point is that most of the 
profiles in the group covered by Figures 3.37 to 3.45 generally show weaker 
Gulf Stream intensity than the moan model of Kcondzhyan. This is probably 
due to the weakened flow in the fall as reported by Isclin (1940) and Fuglister 
(1951). 
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Fig. 3.2 Surface profile based on hydrographic data. Dashed line Keondzayan Model. 





November 30 to December 5, 1932 



Fig. 3.5 Surface pr^'.-le based on hydrographic data. Dashed line Keondzhyan Model. 






August 10-12, 1953 



Fig. 3.7 Surface profile based on hydrographic data. Dashed line Keondzhyan Model. 













Fig. 3.13 Caption see Fig. 3.9 



November 13-15, 1966 



Fig. 3.14 Surface profile based on hydrographic data. Dashed line Keondzhyan Model. 






Fig. 3.16 Caption see Fig. 3.9 




Fig. 3.17 Surface profile bared on hydrographic data. Dashed line Keondzhyan Model. 




October 20 to Novenber 9, 1966 



Fig. 3.20 Surface profile based on hydrographic data. Dashed line Keondzhyan Model 





Fig. 3.22 Caption see Fig. 3.9 




Fig. 3.24 Caption see Fig. 3.9 



Fig. 3.25 Caption see Fig. 3.9 





Fig. 3.26 Caption see Fig. 3.9 






Fig. 3.29 Caption see Fig. 3.9 






Fig. 3.32 Caption see Fig. 3.9 



Fig. 3.33 Caption see Fig. 3.9 





Fig. 3.34 Caption see Fig. 3.9 
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Fig. 3.35 Caption see Fig. 3.9 






Fig. 3.36 Caption see Fig. 3.9 




Fig. 3.37 Caption see Fig. 3.9 
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Fig. 3.40 Caption see Fig. 3.9 




Fig. 3.41 Caption see Fig. 3.9 



Fig. 3.A2 Caption see Fig. 3.9 
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Figure 3.49 Surface profile based on mass transport measurements 
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4. GEOIDAL UNDULATIONS 

4.1 BACKGROUND 

We are interested In obtaining surface elevation information from satellites, 
and any departures from a flat surface over oceans are Interpreted as having 
been caused either by current systems or to a larger extent, by changes in the 
local gravitational acceleration. The continental margins are regions of rather 
dramatic departure of the value of real g from a latltuclinally based constant. 
Reasons for this fact are several, the two foremost due to the changes in ocean 
bottom depth of the order of three to six kilometers over lateral distances of 
the order of tens of kilometers. Factors influencing departures of sea surface 
elevation from a level of constant z are then geostrophic and hydrostatic 
adjustment and local changes in the value of the gravitational acceleration 
coefficient. 

As it occurs, the Gulf Stream between Florida and Cape Hatteras generally 
follows the continental shelf contour. Thus a region of significant changes 
in the local gravitational potential, as is the shelf margin south of Cape 
Hatteras, can also be a regime of Intense sea surface adjustment due to an 
Intense western boundary current. Because of these Integrated complications, 
surface slope calculations cannot be based simply on hydrographic data as done 
in section 2 of this report, but must also be adjusted with geophysical gravity 
determinations. The satellite cannot differentiate between anomaly source, 
it simply sees the resultant. 

Undulations in the geoid and deflections of the vertical can be caused 
by unknown disturbing masses which create gravity anomalies. If one were to 
have complete knowledge of the gravity anomalies , g , around the earth, then 
one could accurately know the other quantities worldwide. There has been a 
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35% coverage, on a 1® by 1® basis, for the entire surface of the Earth. This 
Is a credit to investigators of the twentieth century but unfortunate to the 
determinations of geoidal undulations and deflections of the vertical. Onv‘ 
could of course extrapolate from the existing gravity net to unsurveyed areas 
of the world but confidence in such an extension could not be too high. 

Furthermore, based on our 35%, 1® by 1" coverage, and with the knowledge that 
requirements for anomalies are in area means and point values we can see that 
1® squares may only yield general, as distinct from local, predictions, at 
best. The case of point by point prediction is related to the determination of 
area means of all point anomalies within the area being surveyed. In effect 
then, a large number of points for an area mean or rm actual point observation 
are necessitated when dealing with point predictions. 

Statistical methods have been used to extend gravity data to the unsurveyed 

areas of the Earth. It is generally felt that Ifag is unknown over distances 
3 

between 10 to lO” kilometers then various modifications of statistical methods 
will allow one to fill the gap with reasonably good confidence, given no real 
dramatic change in g across the interval that will not be accounted for in an 
extrapolation. Since, however, we know that gravity anomalies exist, we can 
presume that they are caused by topographic mass deficiencies in the oceans and 
excesses on the continents and by Isostatic compensating surpluses and deficiencies. 
Were the Earth in hydrostatic equilibrium, the equipotential surfaces, including 
level, would be a flattened mathematical ellipsoid (Heiskanen, 1966). It is 
known from geophysical methods though, that the Earth is to appreciable measure 
in. Isostatlc equilibrium (Heiskanen and Meinesz, 1953) so that topographic mass 
excesses and surpluses are compensated for by appropriate masses. Simple response 
of the geold to topographic change would mean that a 5 kilometer high coastal 
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mountain range parallel to a 5 kilometer deep abyssal plain would result in a 
geold undulation of the order of 800 meters. Such changes have not been 
observed either by land, sea or alrbourne Instruments, so there is ho doubt that 
isostatlc compensation does occur to obvious measure. It can be assuredly 
claimed herein that topographic and compensating masses account for not less 
than 85% of g. 

4.2 THEORY AND CALCUUTION 

With U, the potential of the ellipsoid, W, the potential of the geoid, 
r, the gravity on the ellipsoid, g, the gravity on the geoid and ^g, equal 
to g minus t - , the gravity anomaly then the condition that the potentials of 
the ellipsoid and geoid are equal is 

U - W *« Uv + V = constant, (3.1) 

where Uv is the potential of the geoid if there were no mass anomalies and V 
is the disturbing potential. Now, though the potentials of the geoid and 
ellipsoid are equivalent, the two equipotentlal surfaces will not necessarily 
coincide, in which case they will be separated by N, the geoidal undulation in 
units of length. With the assumption that the geoidal undulation N is much 
less than the ellipsoidal geocentric radius, it can be assumed that 

N - ^ (3.2) 

G. G. Stokes developed an expression for V (Heiskanen and Moritz, 1967) 
which when substituted into 3.2 results in the following expression for N 

N ■ — — T J J S(zS') sin-vj/ d“Y doc , (3.3) 


where 
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= BcosN*- fosin^ + 1 -3cos''^( ln( ) J , (3.4) 


and ^ is the angular distance between the area where a g Is determined for to 
the point for which a value for N is being sought, and o( is the azimuth from 
the point causing the effect to the affected point. As is obvious from 3.3, 
integration is carried out over the entire globe, so one needs to have a 
knowledge of the kernal^g, worldwide, in order to perform the integration in 
3.3 The local value of a g is essentially the difference between the locally 
rkhnorved g, g^, minus the international gravity foimuia (Kuiskanen, 1926) 

r- 978.0490 (1+0.005 2884 sin^9 - 0.0000059 sin^ 20 (3.5) 

i 

after several corrections have been applied to the observed gravity readings. 

• It should be noted that the major contributions to N, are functions of values 
of a g within three or four hundred kilometers of the point in question. Also 
of note is the fact that 3.3 is valid under the assumption that all masses 
producing the equipotentlal surface are co.itained inside of the geold. 

The following is a description of marine gravity theory given by Banks 
(1972); and modified by the authors of this report: 

In most areas of the world, Bouguer anomalies are negative in 
mountainous regions and positive over the oceans. These anomalies 
are explainable on the assumption that, above a depth of compensation, 
the average rock density under mountains is less than the average rock 
density under land masses at sea level, while the average rock density 
under oceans is greater than that under land masses at sea level. 

Thus, the effect of topography would be counterbalanced by mass or 
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density differences, manifested as roots or antiroots at the base of the 
earth'? crust. These ideas have led to the hypothesis of Isostatlc 
equilibrium (see e.g., Helskanen and Venlng Melnesz, 1958). 

Although different concepts on the mechanics of Isostac.lc equili- 
brium have been advanced, only two will be used in the calculation of 
isostatic anomalies herein. These are the Pratt-Hayford and 
the Airy-Heiskaiicn concepts of Isostasy, which have been extensively 
used in gravimetric studies. 

The Pratt-IIayford concept assumes that equilibrium is attained 
through uniform vertical columns of different densities, all extending 
downward to the same depth of compensation. Columns in mountainous areas 
are assumed to consist of lighter material than coliuans in the lower 
elevation plains or ocean basins. At the depth of compensation all 
masses are in hydrostatic equilibrium, that Is, they are subjected to 
equal pressures from all directions, whether the masses are under mountains 
or oceans. Thus, two vertical columns of equal cross sectional area, 
extending from the surface of the earth to the depth of compensation, will 
exert the same weight at the base of the columns. 

The pressure at the base of Column 1 equals that at the base of 
Column 2 . 

Then, 

<f“ P')8l (D + “f«2® 

where p is the density of Column 2, which has a surface elevation of sea 
level, f-f' is the density of Column 1, which has a surface elevation 
other than sea level, g^^ is the average value of gravity over Column 1 
and ^2 average value of gravity over Column 2. D is the depth of 
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compensation below sea level and h represents the elevation above sea 
level. » the density Increase or decrease of each column, is a function 
of *'he topographic elevation h, the density , and the depth of compensation 


D. 

T!;e Pratt-Hayford concept also assumes that each column is an Indepen- 
dently compensated unit, that is, columns do not include compensations 
for topography In adjacent columns (see e.g., Helskanen and Vening Melnesz, 
1958) . 


Equation 3.6 shows that at the depth of <:am|-ensatlon, pressures at 
the base of all columns are Identical, and thus the weights of the total 
columns are the same. This Is true assuming that average gravity Is the 
same over both columns which Introduces a small error, as gravity decreases 
with elevation. 

Then 


f 


D 

D+ h 


(3.7) 


which reduces to: 


1 4- 




and, by neglecting higher order terms 

f ^ Jl_ 

f D 


(3.8) 


(3.9) 


gives 



(3.10) 


Isostatlc equilibrium can be expressed mathematically a-3 equality of 
pressure at the depth of compensation or as equality of mass In overlying 
columns. Eqr tlon 3.10 does not satisfy either equality of pressure or of 
mass, but approximates the equality of pressure more closely. This Is due 
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to the fact that only a hwiII change u£ gravity is neglected for the 
different heights of the columns. According to Helslcrmen and Venlng 
Meinesz (1958) "the comcion opinion today is that Isostatlc compensation 
seldom occurs according to the Pratt-Ilayford assumption." 

The Airy-Heiskanen Isontatic concept assumes that a lighter crust 
floats on a heavier mantle, much as an iceberg floits in water. An 
Inadequacy in thl'; (.''n' - ja tha;- r'. .. root of rn Iceberg extends 

downward from sea level, '-jn lo the root oi a i'-wl 'iiass extends downward 
from an average depth of the crust (depth oL sensation), usually 
assumed to be 20, 30, 40, or 60 km below sea level. 

In the Alry-lleiskanen system, materials that compensate for topographic 

variations are assumed to lie vertically beneath the topographic features. 

The density of the entire crust is assumed to be a uniform 2.67 gr/cm^, 

3 3 

and that of the underlyln? mantle 3.27 gr/cm . or 0.6 gr/cm greater than 
that of the crust. The density of 2.67 gr/cm for the earth's crust is an 
arithmetic mean of six estimates that were made between 1811 and 1882 
(Barkness, 1891) this figure has been used by later investigators such as 
Bayf ord and Bowie (1912) . 

Vhen calculating Isostatlc anomalies according to the Airy-Heiskanen 
system, the thickness of the crust for a sea level elevation (land area) , 

T, must be assumed. As thicknesses of the crust are not exactly known, 
Isostatlc anomalies are often computed using various values of T. Analyses 
of anomalies based on these different values of T can provide indications 
of the probable thickness of the crust. 

For Isostatlc anomalies at sea level stations, two corrections are 
applied. The first is the topographic correction, which removes the effect 
of all masses above sea level. This correction is also used to replace 
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ocean water by material having a density equal to the mean surface density 
of the solid earth. There are various methods for obtaining topographic 
corrections, based on reading elevations from topographic maps for certain 
zonal areas (see e.g., Gutenberg, 1959). Although the topographic 
corrections can be made for the entire earth, -'t is sufficient in some 
regions to correct for topography at distances out to 1^6.7 km of the 
station. 

The second correction is made for the effect of the compensating mass, 
which is presumed to be directly beneath each topographic feature and 
below the assumed depth of compensaclor Calculation of this second 
correction Is lengthy; published tables and maps (Hayford and Bowie, 1912; 
Helskanen, 1938; and Karkl, Kivioja, and Heis’:anen, 1961). 

The Bouguer anomaly for stations at sea level Includes only a correction 
for the effect of the surrounding topography. The topographic correction 
used In the Isostatic anomaly calculation is also the Bouguer topographic 
correction. These topographic corrections are always added to thu observed 
sea level gravity observation. This is due to the fact chat the mountains 
above the station and the density deficiency of the water below the station 
both have the effect of reducing the observed gravity measurement from 
what would have been measured at sea level with no water or topography 
present. 

The Bouguer reduction is sometimes viewed as a type of Isostatic 
reduction corresponding to an infinitely thick crust (Helskanen and 
Venlng Heines z, 1958). 

Positive or negative isostatic anomalies are Interpreted to mean that 
hydrostatic equilibrium does not exist at the assumed depth of compensation. 


I 



or thac the method of obtaining the anomalies needs some modification 
(Coulomb and Jobert, 1963). 

Isostatlc theories propose that at some depth below sea level the 
pressure on equal areas is the same everywhere. The horizontal dimensions 
required for a region to be in equilibrium are not defined nor generally 
agreed upon. Some Investigators believe that cross sections of the 
order of 10,000 kra may be large enough to be compensated (Bowie, 1917); 
other studies indicate that an area the size of the Hawaiian Islands 
(16,638 km2) is too small to be In complete equilibrium (Niskanen, 1945). 
Some people envision that an area must be as large as the United States 
before equilibrium can be achieved (Bowie, 1924). Tsuboi (1940) has 
calculated that equilibrium requires an area with the horizontal dimension 
about three times the depth of compensation. This estimate is generally 
accepted as reasonably correct. Pu** to the strength of the earth's crust, 
it is believed that complete equilibrium is extremely unlikely for any 
topographic feature (Heiskanen and Vening Meinesz, 1958). 

The depth of compensation or the proper thickness of the crust must 
also be assumed when calculating isostatic anomalies. Much as Bouguer 
anomalies can be considered to be as forms of isostatic anomalies when 
the depth of compensation is infinite (completely rigid materials of 
great thickness) , so free-air anomalies can be considered as forms of 
Isostatic anomalies when the depth of compensation is at sea level (zero 
crustal thickness, with topographic masses having zero rigidity). A 
geologically more realistic picture, in view of the high but not infinite 
rigidity of crustal rocks, is that there is a depth of compensation at 
some real depth beneath the surface. The depth at which compensation 
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occurs Is assumed. In the Airy-Hclskanen concept, the assumed depths 
are 20, 30, AO and 60 km; in the Pratt-Hayford concept, the assumed 
depth is 113.7 km. 

A particular density structure for the crust and mantle is also 
assumed for the various isostatic systems. Although the constant 
crustal density assumed in the Airy-Heiskanen theory may closely approxi- 
mate the world averaj^e, specific areas of the world may deviate from 
this assumed value (see e.g., Walcott, 1967). This fact could affect the 
results obtained from an Isostatic study. 

In making Airy-Heiskanen isostatic reductions, the Mohoroviclc 
discontinuity is assumed to be the boundary between the crust consisting 
of one uniform density and the mantle consisting of another uniform 
density. Seismologists find that two and possibly three discontinuities 
exist above the Mohoroviclc discontinuity, and possibly some beneath the 
discontinuity. Calculated Isostatic anomalies have shown, however, that 
nearly the same Isostatic result can be obtained by assuming either one 
or even several discontinuities. Thus, before two or more discontinuities 
can be used in calculations, more reliable evidence is needed concerning 
the discontinuities (Heiskanen and Vening Meinesz, 1958). 

So, one can correct gravity observed at sea by making four corrections 
(Grant and West, 1965). The first correction is the Free-Air correction, 
which takes into account the position of the gravimeter with respect to sea 
level. The correction is subtracted from g^ as a function of depth of measure- 
ment below sea level, and is given by the following equation 

Agj^ - 0.3086d. (3.11) 
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This correction, of Itself at sea, has little meaning because of the fact of 
having to account for the thin plate of water above the ocean floor. The Bouguer 
Plate Correction must be applied herein. This correction deletes the effect 
of the upward pull of the water section and in a second application will restore 
the pull downward by the water section if the measuring device had been located 
at the ocean surface. This reduction Is always positive in nature and is given 
as 

Agj^^ - O.OSeid. (3.12) 

It should be mentioned that the combined corrections to g^, as given by Ag^^ 
and Agj^ are referred to as the free-alr-at-sea correction. 

A correction, based on the assumption that the ocean basins represent 

a mass deficiency not compensated isostatically , and is then always positive. 

This correction is called the Secondary Plate Correction and is given by the 
following relation 

(3.13) 

where h Is total ocean depth. 

A fourth correction, Ag^, called the Isostatic Correction can be applied 
and can be based on the aforementioned theories of Pratt-Hayford and Airy- 
Helskanen. 

Another correction is that due to the fact that measurements are being made 
from a moving vehicle and is called the Eotvos correction and must be removed 
from observations. It is 

^®ec " * 0.14584 Vg^, (3.14) 

where the plus sign denotes eastward ship movement, the negative sign denotes 
westward motion and Vg^ denotes the ships velocity. This correction is based 
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on the fact of existence of the ficticious centrifugal force due f) tin Earth's 
rotation. 

Knowing the constituents of the gravity corrections then, oae derives a 
series of resultant anomalies, representing the free-air-anomaly-at-^ea, Ag^ 


the Bouguer anomaly at 

sea 

, Ag^ and the isostatlc 

anomaly at sea, Ag^. 

The 

appropriate relations 

are ; 

as follows 



Agf 

»o 

" ^Sfc ^'^bc ■ ^^ec ~ 


(3.15) 



“ Agf + 


(3.16) 

and 

ASi 

= Agj^ + Ag^^ 


(3.17) 


From comparisons, it appears that the free-air anomalies, are generally 

like those which result from isostasy reduction considerations and c'e then, 
considering the greater simplicity of calculating g^, used to calculate the 
Stokes function N in expression 3.3. 

Worzel and Shurbet (1955) compiled a series of gravity observations made 
at 104 sea stations. Three of their seven transects (Figure 4.1) are appropriate 
to our Gulf Stream area of interest and the marine gravity data which has been 
used in an attempt to calculate the Stokes function is shown in Figures 4.2 
through 4.5. 

It has been mentioned previously that, as Garland (1965) said, the calcu- 
lation of N at any point is dependent on a knowledge of Ag over the whole earth 
and there are barely sufficient observations in some regions to permit accurate 
calculations anywhere. To avoid this problem, then, Vening Meinesz (Heiskanen 
and Meinesz, 1958) demonstrated that where gravity anomaly gradients are "smooth," 
then the azimuth of the near field vertical deflection is measurably co-incident 


i 




with the direction of the gravity anomaly gradient. If one draws a circle 
of radius 30 kilometers (e.g., Heiskanen and Heines?., 1958) around the point 
In question and computes average free-air gravity anomaly giadiencs across the 
meridional and zonal diameters of the circle then, one can use the Mcinesz 
development to imply an approximated Stokes function. This proced ‘a was 
followed in the development of the curves depicted in Figure 4.6 th,-ough 4.8 
and was based on the hctzcl and Shurbet d.'.ta. 

The results, as shown in the curves indicate that in general the Stokes 
function Increases in magnitude from Cape May to Cape Henry to Cape Hatteras 
with the total differential in the geoldal undulation being augmented from 
approximately one meter at Cape May to two and a half meters at Cape Henry to 
three meters at Cape Hatteras. 

Since the order of magnitude of these undulations is at least as large 
auo usuaxly several times as large aa sea level cha*'gc due to Culf Stream 
systems adjustment (cf. Sections 2 and 3 of this report) then the non-ocean 
vturrent based causes must be subtracted from the total deviation in the geoid 
if one is to look at sea surface slopes related to ocean systems with any 


Insight and confidence. 
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5. DISCUSSIONS AND RECOMMKNDATIONS 

5.1 DISCUSSIONS 

Based on the results of this study, the foMowlng points are of special 

Interest: 

(a) Contrary to the accepted view, there ate fluctuations of the Gulf Stream 
occurring over extremely short periods of time of the order o^’ several 
days. Such changes can easily be observed from Figures 3.32 to 3.36. 
Qualitatively, such changes are independently confirmed by infrared 
imagery as shova by Stumpf (197A) in Figures 5.1 and 5.2. It is not clear 
what the cause of such short time change may be. Local wind stress seems 
improbable both from the energy argument of Stomcal (1966) and from actual 
field data presented in Chapter 3 which shows iow correlation between 
local wind and the Gulf Scream location and strength. A possible cause is 
the Instability of the stream which causes some Guit Stream water masses cu 
spin off in an apparent energy relieving process. Such "eddies," "shingles" 
or "sausages" would show up in hydrographic data as well as infrared imagery 
ns reported by Rao, Strong and Koffler (1971). These phenomena are 
impossible to study using only shipboard instruments and methods. This 
can be further demonstrated by the following example. Lat curve 1 in 
Figure 5.3 be the Gulf Stream at time Tj^, and curve 2 be the Gulf Stream 
at time T 2 * If a ship set cut to study the phenomenon at Tj and finished 
the study at T 2 . The data that the ship actually collected would result in the 
solid curve which is not true at any time. Of course if the time scale of 
the phenomenon Is relatively long compared to the cruise time, the ship's 
data will be more accurate. Unfortunately, high frequency changes do 
occur. Such mobility of the Gulf Stream immediately leads to the question 
not only of the fidelity of the data, but also of the geostrophic assumption 
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which depends crucially on steady state conditions. The Influence of 
short period variations on dynamic calculations has been discussed by 
Seiwell (1939). The implication seems to be that for highly variable, 
large scale motions, such as world current systems, remote sensing may be 
the only practical tool to use in a synoptic sense. 

(b) Tlie seasonal variations of the Gulf Stream are considerable. Due to the 

seasonal changes of t e global wind field, the Gulf Stream ill have to 
change accordingly. This problem has been studied by Fuglister (1951) and 
he concluded: "The maximum current speeds occur during the summer in the 

southern segments of the Gulf Stream system and during the winter in the 
north segments. The minimum speeds occur during the fall months throughout 
the system." The results are presented graphically in Figures 5.4 and 
5.5. The fluctuations are, of course, due to the global wind field and 

are not directly related to the local wind field. Fuglister found, however, 
that the Gulf Stream strength correlated directly with the trade winds but 
only partially with the local mean wind. 

(c) Beyond the seasonal fluctuations, nothing can be said about year to year 
changes. The existing data are too sparse to be of any use in addressing 
long term variational patterns. However, since long .term meteorological 
cycles do exist, it stands to reason that there may be similar cycles 
which the Gulf Stream observes. Fluctuations at these scales can not be 
studied by traditional methods because for such long term meteorological 
cycles, the temporal scales are of sufficiently large enough size that the 
Global circulation systems (Figure 5.6) might interact with each other. 

Such interaction will, in turn. Influence long term weather cycles. This 
phenomenon should be one of the most Important and unique tasks for remote 
sensing oceanography to address Itself to. 
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(d) Due to the discussion above, It becomes obvious that a single Gulf Stream 
model, good for all time. Is Impossible. Though a mean Is still possible, 
as reported by Keondzhyan (1973), the model should be used judiciously. 

(e) The undulations caused by topography and Isostasy In the region of the 
continental shelf break would appear to be of the order of several meters 
excursion which Is at least as large In magnitude as any change In sea 
surface height as a function of the Gulf Stream. Since the solid earth 
based causes for geoldal undulations are not time dependent, then it would 
be of considerable worth to be able to simply delete such cause and effect 
from the data. The Implication is that subtracting out the effects of 
crossing the shelf rise would allow one to look at ocean systems directly. 

5.2 RECOMMENDATIONS 

Based on the study, the following recotmendations can be made: 

(a) The remote sensing method should be used to cover the global ocean continuously 
for an extensive period of time so that the temporal and spatial variations 
can be determined over all ranges. It is granted that there are limitations 

of remote sensing methods and, therefore, comparisons with available in situ 
field data should be made complementarily with each other. 

(b) Since most strong current systems in the world occur as western boundary 
flows, e.g. , the Gulf Stream and Kuroshio, care should be taken when we 
are trying to infer current conditions from altimeter measurements and 
geostrophic theory, since part of the Gulf Stream generally follows the 
continental shelf where the geold perturbation caused by the continental 
rise could be mistaken as part of the stream. Therefore, a detailed study 
of geold perturbation along the continental shelf should be one of the 
most crucial tasks to be worked on. 
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Figure 5.1 Facsimile chart for October 23-28, 1973 (Stumpf, 1974). 
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Figure 5.3 The result of a slow ship survey In a fast changing field 
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Figure 5.4 Location ot segments or the Gult oLi.C(iul SyStCvu 
used to calculate annual variations in current 
speeds. Arrows show resultant current directions 
(Fuglister, 1951). 


I 




80 


70 



lAl 


JUN ri8 APR MAY JON A'L AUS StR OtI t.IV OtC 



FLORIDA 


sovm* 

OF 

HATTERAt 


H.E. OF 
NATTCRAS 


Figure 5.5 Annual variation in the curient 
speed in three segments of the 
Gulf Stream system. Observed 
data, solid lines calculated 
combination of annual and semi- 
annual periods, dashed lines 
(Fuglister, 1951). 
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